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Abstract.  The  results  of  nighttime  measurements  of  positive  ion  densities  between  240  and 
1S75  km  are  presented.  The  ion  scale  height  is  found  to  change  abruptly  around  600  km,  in¬ 
creasing  from  120  km  at  500  km  to  1400  at  640  km;  it  increases  gradually  with  height  above 
this  level.  An  ion  temperature  of  1050  ±  50°K  is  deduced  from  the  ion  scale  height  measured 
above  the  peak  of  the  F  region.  The  ion-density  measurements  were  obtained  by  means  of 
a  spherical  electrostatic  analyzer  flown  on  a  Blue  Scout  rocket.  In  the  altitude  region  2S0  to 
350  km  the  ion-density  values  agree  to  within  15  per  cent  with  the  electron-density  measure¬ 
ments  obtained  from  local  ground-based  ionosondes.  The  vertical  distributions  of  the  principal 
light  ions  above  400  km  are  calculated,  assuming  the  existence  of  diffusive  equilibrium.  The 
relative  abundance  of  the  different  ionic  species  at  the  reference  level  is  determined  from  the 
measured  profile  of  the  total  ion  density.  The  results  of  the  nighttime  flight  show  that  the 
atomic  oxygen  ion  is  the  principal  constituent  below  530  km  and  ionized  hydrogen  the  major 
constituent  above  550  km.  An  upper  limit  is  placed  on  the  helium  ion  distribution,  which 
attains  its  maximum  value  at  550  km.  At  this  level  it  constitutes  a  maximum  of  about  7  per 
cent  of  the  total  ion  density.  The  ion-atom  interchange  coefficient  for  the  production  of  NO* 
is  calculated  from  the  variations  with  time  of  charge  densities  in  the  altitude  region  280  to 


300  km ;  the  coefficient  is  found  to  vary  between 
Introduction 

By  the  method  described  in  an  earlier  pa¬ 
per  [ Sagalyn  et  al.,  1963]  positive  ion  densi¬ 
ties  have  been  measured  over  the  altitude  region 
240  to  1860  km  during  a  nighttime  flight  of  a 
Blue  Scout  rocket.  A  third  mode  of  operation 
was  added  to  obtain  new  information  about  the 
positive  ion  energy  distribution.  The  results  are 
compared  with  day  and  nighttime  electron-den¬ 
sity  profiles  obtained  by  other  techniques.  Ion 
temperature  and  the  average  ion  mass  and  scale 
height  are  determined  as  a  function  of  altitude 
from  the  experimental  results. 

Nicolet  [1961], Hanson  [1962], Bower  [1962], 
Bourdeau  and  Bauer  [1963],  Ulwick  and  Pfister 
[1963],  Mange  [1960,  1961],  Donley  [1963], 
and  others  have  presented  theoretical  discus¬ 
sions  or  experimental  results  on  charged  particle 
distributions  above  the  ionization  maximum  in 
the  F  region.  Applying  equations  derived  by 
Mange  [1960],  Hanson  [1962],  and  Bauer 
[1962]  we  have  calculated  the  altitude  distribu¬ 
tion  of  the  individual  light  ion  species  above  the 
level  where  diffusive  equilibrum  may  be  assumed 
to  exist.  The  reference-level  composition  is  de¬ 
duced  from  the  analysis  of  the  vertical  distribu¬ 
tion  of  the  total  ion  concentration  obtained  on 
the  Blue  Scout  flight.  The  vertical  profiles  of 


2  and  4  X  10*1J  cm3  sec*1. 

atomic  oxygen  and  hydrogen  ions  are  obtained 
for  altitudes  above  400  km;  an  upper  limit  is 
given  for  the  helium-ion  distribution.  The  re¬ 
sults  are  compared  with  other  nighttime  meas¬ 
urements. 

Langstroth  and  Hasted  [1962],  Talrose  et  al. 
[1962],  Galli  et  al.  [1963],  and  Sayers  and  Smith 
[1963]  have  carried  out  laboratory  measure¬ 
ments  of  the  ion-atom  interchange  rate  coeffi¬ 
cient  for  the  formation  of  NO+,  an  important 
ion  loss  process  in  the  nighttime  F  region.  The 
applicability  of  their  results  to  the  ionosphere 
is  not  certain,  as  the  experiments  usually  pro¬ 
vide  only  an  upper  limit  for  the  reaction;  also, 
the  reactants  in  the  laboratory  and  ionosphere 
may  not  be  in  the  same  states.  Yonezawa  et  al. 
[1959],  Bates  and  Nicolet  [1960],  Danilov 
[1960],  and  others  have  calculated  values  for 
this  reaction  from  upper  atmospheric  data 
under  various  assumptions.  The  results  vary 
by  three  orders  of  magnitude. 

The  ion-atom  interchange  coefficient  for  the 
formation  of  NO+  is  calculated  in  this  report 
using  nighttime  variations  of  charged  densities 
obtained  from  the  Blue  Scout  data  and  electron- 
density  profiles  from  a  local  ionosonde  station. 
The  continuity  equation  is  solved  under  the 
assumption  that  between  280  and  300  km  the 
time  variation  of  positive  ions  is  principally  due 
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Fig.  1.  Simplified  block  diagram  of  spherical  elec¬ 
trostatic  analyzer. 


to  ion-atom  interchange,  diffusion,  and  radiative 
recombination. 


The  Experiment 

A  simplified  block  diagram  of  the  spherical 
electrostatic  analyzer  experiment  is  shown  in 
Figure  1.  This  experiment  is  an  extension  of 
the  spherical  ion  probe  first  used  by  Griiigauz 
and  Zelikman  [1957]  and  Gringauz  et  al.  [1961] 
for  satellite  measurements.  The  experiment  is 
operated  in  three  modes. 

Mode  1.  A  fixed  saturation  voltage  is  ap¬ 
plied  between  grid  1  and  the  collector.  As  shown 
in  Sagalyn  et  al.  [1963]  the  ion  density  N  is 
related  to  the  collector  current  i  through  the 
relation 


i  =  ^  r2  Necj(  V)a 

A 

•  [>  +  irU2{x  +  erf  (x) 

where 


(1) 


r  —  radius  of  grid  1. 
e  =  ion  charge. 

c  =  most  probable  ion  velocity. 
q  =  rocket  velocity. 

a  =  experimentally  determined  transmission 
factor  for  grid  1. 
x  =  2  q/cT11*. 

/(V),  a  function  of  rocket  potential  V,  is 
obtained  from  mode  2  analysis. 


2  fz 

erf  (x)  =  -175  I  exp  ( —  if)  dy. 

7T  t/0 


Mode  2.'  Saturation  voltage  is  maintained 
between  grid  1  and  the  collector  while  a  linear 
sweep  voltage,  range  ±  17  volts,  duration  3  sec¬ 
onds,  is  applied  to  grid  1.  The  rocket  potential 
Vr  and  }(V)  of  equation  1  are  obtained  from 
the  resulting  current  voltage  curves.  Current 
voltage  curves  obtained  at  436  and  1650  km  on 
a  Blue  Scout  flight  launched  on  April  12,  1961, 
are  shown  in  Figure  2.  The  analysis  of  this 
mode  is  given  in  Sagalyn  et  al.  [1963],  where 
it  is  shown  that 

VR  =  (m/2e)(q  +  2.3 c)2  -  AF,  (2) 
where 

VR  is  the  rocket  potential  with  respect  to  the 
undisturbed  plasma. 
m  =  mass  of  heaviest  ion  in  the  gas. 
c  =  most  probable  ion  velocity. 

AFi  =  value  of  sweep  voltage  corresponding  to 
complete  retardation  of  most  energetic 
ions  in  the  gas. 

Values  of  f(V)  and  VK  measured  as  a  function 
of  altitude  on  the  Blue  Scout  flight  are  given  in 
Table  1. 

The  analysis  of  mode  2  data  to  obtain  ion 
mass  and  temperature  will  be  presented  in  a 
future  paper. 

Mode  3.  On  the  flight  under  discussion  the 
spherical  analyzer  was  operated  in  a  third  mode 
(see  Figure  1),  in  which  the  voltage  between 
grid  1  and  the  collector  was  varied  in  steps  be¬ 
tween  0  and  100  volts.  From  the  resulting  cur¬ 
rent  voltage  curves  the  assumptions  made  in 
the  calculations  of  the  required  saturation  volt¬ 
age  [Sagalyn  et  al.,  1963]  can  be  tested.  The 
voltage  that  must  be  applied  between  the  elec¬ 
trodes  for  collection  of  all  positive  ions  in  the 
gas  may  be  written 

V  =  (E/em’/a2)  -  1]  (3) 

where  E  =  energy  of  the  most  energetic  ion  in 
the  gas  striking  grid  1  at  grazing  incidence,  b  — 
radius  of  grid  1,  a  =  radius  of  the  collector. 

Mode  3  current  voltage  curves  obtained  at 
various  altitudes  on  the  Blue  Scout  flight  are 
shown  in  Figure  3.  At  617  km  the  saturation 
potential  was  measured  to  be  55  volts;  using 
equation  3,  and  assuming  a  composition  of  60 
per  cent  0+,  40  per  cent  Il+,  it  is  calculated  to 
be  50  volts.  As  shown  in  Figure  3,  the  satura¬ 
tion  voltage  was  found  equal  to  33  volts  at  1800 
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Fig.  2.  Collector  current  versus  sweep  voltage  applied  to  grid  1;  mode  2;  altitudes  436  and 

1650  ivm. 


km;  it  was  calculated  to  be  28  volts  assuming 
II+  the  dominant  ion. 

A  differential  energy  distribution  for  charged 
particles  in  the  gas  is  also  obtained  from  these 
data. 

Results  and  Discussion 

(a)  Vertical  Distribution 

The  vertical  distribution  of  positive  ions  meas¬ 
ured  on  the  ascent  and  descent  portions  of  a 
Blue  Scout  flight  launched  from  Cape  Canaveral 
on  April  12,  1961,  at  0107  EST  is  shown  in  Fig¬ 
ures  4a  and  46.  The  maximum  ion  density  is 
measured  to  be  2.25  X  10!  ions/cm*  at  340  km 
on  the  descent.  This  is  in  good  agreement  with 
the  ionosonde  value  of  2.4  X  10"  electrons/cm' 
measured  at  the  nearest  ionosonde  station  at 


San  Salvador.  At  300  km  on  the  descent,  the 
latitude  of  the  Blue  Scout  rocket  is  2°  south 
of  the  San  Salvador  station.  The  ionosonde  data 


TABLE  1.  Altitude  Variation  of  Rocket  Potential 
Vr  and  f(V) 

Altitude,  V„ 

km  /(V)  volts 


300  0.98  ±0.2  -1.4  ±0.3 

500  1.0  -1.3 

700  1.0  -1.5 

900  1.1  -1.4 

1100  1.2  -1.5 

1300  1.3  -1.8 

1500  1.5  -2.0 

1700  1.5  -2.2 

1875  1.6  -2.1 
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ION  DENSITY  (No./cm> 


Fig.  4a.  Altitude  versus  positive-ion  density,  Blue  Scout  flight,  0102  EST,  April  12,  1961, 

Latitude  28°  N,  descent. 
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Fig.  4b.  Altitude  versus  positive-ion  density,  Blue  Scout  ascent;  ionosonde  electron-density 

data  obtained  at  latitude  28°N. 


shown  in  Figure  4a  were  obtained  from  G.  Stone- 
hocker,  of  the  National  Bureau  of  Standards, 
Boulder,  Colorado.  A  sharp  change  in  the  slope 
of  the  distribution  is  observed  near  700  km  on 
the  ascent  and  560  km  on  the  descent  portion  of 
the  flight.  Above  600  km  the  descent  values  are 
about  30  per  cent  higher  than  the  ascent  values; 
below  this  level  the  ascent  and  descent  values 
agree  within  10  per  cent.  The  difference  between 
ascent  and  descent  values  should  be  small  in  the 
F  region,  since  the  total  change  in  geomagnetic 
latitude  for  the  /'’-region  data  is  only  4°. 

Above  700  km  the  vertical  profiles  differ  sig¬ 
nificantly  from  other  profiles  recently  presented 
[Hanson,  1962;  Burning,  1960;  Gringauz,  1962], 
This  will  be  discussed  further  in  the  next  section. 
The  results  and  discussion  in  (a),  (6),  and  (c) 
were  presented  orally  by  Sagalyn  et  al.  [1962] . 


(6)  Ion  Temperature,  Average  Ion  Mass 

Assuming  diffusive  equilibrium  to  exist  above 
300  km,  the  ion  temperature  is  calculated  from 
the  hydrostatic  equation  for  charged  particles 
(equation  4)  to  be  1025  ±  50°K  on  the  descent 
and  1100  ±  50°K  on  the  ascent.  Mange  [1960] 
derived  equation  4  assuming  that  the  electron 
and  ion  temperatures  arc  equal.  The  experi¬ 
mental  observations  of  Evans  [1962]  indicate 
that  this  is  a  valid  assumption  for  the  nighttime 
F  region.  Brace  et  al.  [1963]  recently  reported 
results  of  a  nighttime  flight  in  which  the  ascent 
and  descent  electron  temperatures  differed  by 
5  to  15  per  cent  in  the  altitude  region  280  to  350 
km.  They  indicated  that  this  result  implied  that 
the  electron  and  ion  temperatures  were  not 
equal.  Assuming  their  interpretation  to  be  cor¬ 
rect,  the  effect  on  the  numerical  solution  of  the 
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hydrostatic  equation  is  negligible.  On  the  basis 
of  the  measurements  of  Istomin  [1960]  it  was 
assumed  that  0+  was  the  dominant  ion  between 
300  and  500  km. 

Assuming  the  atmosphere  to  be  isothermal 
above  500  km,  the  variation  with  altitude  of  the 
average  ion  mass  and  scale  height  can  then  be 
calculated  [Mange,  1960] 


-log.  n(z)  — 


gm+(z) 

2kT 


(4) 


where 

g  —  acceleration  of  gravity  =  goR1/ (B  +  Z)'\ 
g„  is  the  acceleration  at  R,  the  earth’s 
surface. 

k  =  Boltzmann’s  constant. 

T  =  ion  temperature. 

m+{z)  =  average  ion  mass. 

Hi  =  ion  scale  height. 

n{z)  =  total  ion  density  at  altitude  Z. 

The  results  obtained  from  the  descent  data  are 
shown  in  Table  2.  A  rapid  decrease  in  the  aver¬ 
age  mass  is  found  between  600  and  800  km; 
this  is  qualitatively  indicated  in  the  vertical 
profiles  of  Figure  4  showing  a  rapid  change  in 
slope  in  this  region. 

The  diffusive  equilibrium  equation  must  be 
used  with  caution  in  the  interpretation  of  day¬ 
time  measurements,  since  it  has  been  shown  by 
the  experimental  work  of  Spencer  et  al,  [1962] 
and  Brace  et  al.  [1963],  and  by  the  theoretical, 
work  of  Hanson  [1963]  and  Dalgarno  [1963], 
that  large  differences  can  exist  between  the  elec¬ 
trons  and  the  gas  temperatures  in  the  lower  F 
region  during  the  day.  These  references  also 
show,  however,  that  the  daytime  electron  and 
gas  temperatures  approach  equality  at  about 
330  km  during  quiet  ionospheric  conditions.  The 
diffusive  equilibrium  equation  should  therefore 
be  applicable  above  330  km  under  quiet  daytime 
conditions. 


(c)  Vertical  Distribution  oj  the  Ionic 
Constituents 

Dungey  [1955]  first  discussed  the  distribution 
of  charged  particles  in  the  earth’s  outer  atmos¬ 
phere  under  diffusive  equilibrium  conditions. 
Mange  [1960,  1961],  Bauer  [1962,  1963],  and 
Hanson  [1962]  have  shown  that  if  diffusive  equi¬ 
librium  exists  information  can  be  obtained 
about  the  distribution  of  the  separate  ionic  con¬ 


TABLE  2.  Variation  of  Average  Ion  Mass  m+(z 
and  Ion  Scale  Height  with  Altitude,  Descent 


Altitude, 

km 

amu 

IT.  * 

11  \  > 

km 

400 

16 

1.3  X  102 

500 

17 

1.2  X  1Q‘ 

640 

1.5 

1.4  X  103 

760 

1.3 

2.0  X  103 

920 

1.0 

2.0  X  103 

1000 

1.2 

2.0  X  103 

1200 

1.2 

2.1  X  103 

1400 

1.1 

2.4  X  10s 

1600 

1.0 

2.8  X  103 

1875 

1.1 

2.7  X  103 

*  The  accuracy  of  the  values  in  this  column  is 
15  per  cent. 


stituents  from  the  measurement  of  the  vertical 
profile  of  the  total  charge  density.  Bates  and 
Patterson  [1961,  1962]  calculated  the  distribu¬ 
tion  of  ionic  species  for  model  atmospheres  con¬ 
taining  two  or  three  ionic  species.  Considerable 
interest  in  the  experimental  determination  of 
ion  composition  above  the  F2  peak  has  resulted 
from  Nicolet’s  [1961]  calculations  which  indi¬ 
cated  that  helium  ions  might  be  present  in  sig¬ 
nificant  quantities  in  the  earth’s  upper  atmos¬ 
phere. 

Bauer  [1962]  solved  the  hydrostatic  equation 
to  obtain  the  altitude  distribution  for  the  total 
number  of  charged  particles  of  a  given  sign  in 
terms  of  geopotential  altitude.  Hanson  [1962] 
integrated  the  equation  for  the  ionic  constitu¬ 
ents  and  calculated  the  distribution  of  the  indi¬ 
vidual  ionic  species,  making  estimates  of  the 
ion  composition  at  a  chosen  reference  level. 
Comparing  his  calculations  with  Hale’s  ion- 
density  measurements,  Hanson  concluded  that  a 
helium  layer  approximately  2000  km  thick  must 
exist  in  the  altitude  range  1000  to  3000  km  and 
that  diurnal  variations  of  this  layer  should  be 
small.  Bourdeau  et  al.  [1962]  concluded  from 
the  retarding  potential  analysis  of  a  planar  ion 
probe  at  1650  km  that  the  ratio  He  VO*  was  1.5 
for  a  daytime  measurement. 

In  the  analysis  under  discussion,  the  distribu¬ 
tion  of  the  ionic  constituents  above  the  peak  of 
the  F  layer  is  obtained  without  making  assump¬ 
tions  about  a  reference-level  composition.  As¬ 
suming  diffusive  equilibrium  to  exist  for  charged 
particles,  the  general  expression  for  the  concen- 
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tration  of  charged  particles  of  a  given  sign  n{z) 
is  written  in  terms  of  altitude,  ion  mass,  and 
temperature.  The  relative  proportions  of  the 
ionic  constituents  at  a  reference  level  are  deter¬ 
mined  by  comparing  measured  charge  density 
distributions  with  solutions  of  the  general  ex¬ 
pression  for  n(z)  for  a  wide  range  of  reference- 
level  compositions  and  temperatures.  The  com¬ 
position  and  temperatures  that  are  in  best 
agreement  with  the  experimentally  determined 
profiles  are  used  to  solve  exactly  the  hydrostatic 
equation  for  the  individual  ionic  species.  The 
procedure  is  applied  to  the  results  of  the  Blue 
Scout  flight  to  obtain  the  altitude  variations  of 
0*,  IT,  and  He+. 

Assuming  that  the  gas  is  electrically  neutral, 
the  general  expression  for  the  variation  of  total 
ion  density  with  altitude  is  obtained  from  the 
integral  form  of  the  hydrostatic  equation  for 
charged  particles.  Bauer  [1962]  derived  this  ex¬ 
pression  in  terms  of  geopotential  altitude.  It 
can  also  be  obtained  by  summing  over  the  ex¬ 
pression  of  the  individual  ionic  species  derived 
by  Hanson  [1962] 


n+0)  =  M+(z0)|Z)  «■,■(«„) 

R\z  —  z»)  , 

(R  +  z)(R  +  Zo)_  / 


exp 


'  kf 


12  «i(* o) 


(5) 


where  n+(z0)  is  the  total  concentration  of  positive 
ions  at  the  reference  altitude  z0,  wi<  is  the  mass 
of  the  ith  ionic  constituent,  n,(z0)  is  the  con¬ 
centration  of  the  ith  constituent  at  the  reference 
level  Zo.  Equation  5  can  be  evaluated  without 
knowledge  of  the  individual  ion  densities  rqfzo) 
at  the  reference  level,  since 

n,(z0)  =  a,(z0)  J2  ni(z o)  (6) 

J 

where  at(z0)  is  the  relative  abundance  of  the  ith 
constituent  at  the  reference  height. 

Substituting  n,(z„)  given  by  equation  6  in  the 
numerator  of  equation  5,  the  general  expression 
for  the  charge  density  becomes 


n+(z)  =  n+(z0)  { 12  oci(z0)  exp  —  [i/om.-B2 

i 

•(z  -  Zo)/kT(R  +  z)(R  +  z„)]}1/2  (50 

The  dependence  of  this  solution  on  tempera¬ 
ture  and  ion  composition  is  illustrated  in  Figures 


5  and  6.  In  Figure  5  the  solution  of  equation  5' 
is  given  for  temperatures  ranging  between  800 
and  1600°K  and  for  two  reference-level  compo¬ 
sitions:  0.5  per  cent  IT,  0.5  per  cent  Iie+,  99 
per  cent  0+;  and  5.0  per  cent  H+,  0.0  per  cent 
TIe+,  95  per  cent  0+.  It  is  seen  that  the  lower  the 
temperature  the  more  rapid  the  decrease  of  ion 
density  with  height.  The  magnitude  of  the  de¬ 
crease  in  charge  density  for  a  given  temperature 
change  varies  considerably  with  the  reference- 
level  composition. 

The  dependence  of  the  solution  of  equation  5' 
on  composition  is  shown  in  Figure  6.  The  tem¬ 
perature  is  fixed  at  1400°K  while  the  propor¬ 
tions  of  the  hydrogen,  helipm,  and  oxygen  ion 
are  varied.  As  is  implied  by  equation  5',  the 
greater  the  average  ion  mass  at  the  reference 
level  the  more  rapid  the  decrease  in  total  ion 
density  with  height. 

It  is  seen,  therefore,  that,  if  the  temperature 
is  known  or  can  be  independently  determined, 
bounds  can  be  placed  on  the  composition  at  the 
reference  level  by  comparing  theoretical  and  ex¬ 
perimentally  determined  distributions.  The  re¬ 
sults  of  the  Blue  Scout  flight  are  shown  in  Fig¬ 
ure  7  in  terms  of  z  —  z0  with  z0  =  425  km  for 
the  ascent  and  350  km  for  the  descent  data.  N0 
is  the  total  ion  density  at  the  reference  level  z0. 
The  solution  of  equation  5'  for  the  previously 
determined  temperature  of  1050°K — see  section 
( b ) — and  varying  composition  at  the  reference 
level  is  shown  in  the  same  figure.  It  is  seen  that 
a  reference-level  composition  of  5  per  cent  H+, 
93.5  per  cent  0+,  1.5  per  cent  He*  best  fits  the 
data.  By  means  of  perturbation  techniques,  a 
refined  determination  of  the  reference-level  com¬ 
position  can  be  obtained  taking  the  above  result 
as  first  best  guess.  The  analysis  gives  5  ±  2  per 
cent  H+,  94  ±  2  per  cent  O*,  and  He+  between 
0  and  1.5  per  cent;  1.5  per  cent  represents  an 
upper  bound  for  the  helium  content  at  the  ref¬ 
erence  level.  Since  the  total  ion  concentration  at 
the  reference  level  was  measured  to  be  1.35  X 
108  cm-*,  H+  =  6.75  X  10*  cm-*  and  0*  =  1.27 
X  108  cm"*.  The  solution  gives  a  value  of 
He+  <  1  X  10*  cm"*  at  z0. 

Ionic  species.  Hanson  [1962]  derived  an  ex¬ 
pression  for  the  ith  constituent  in  terms  of  the 
reduced  altitude  z'.  For  purposes  of  computa¬ 
tion  it  is  more  convenient  to  write  this  equation 
in  terms  of  mass,  temperature,  and  altitude, 
above  the  surface  of  the  earth: 


-Zo  (km) 
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Fig.  5.  Dependence  of  ion-density  distribution  on  temperature  for  reference-level  composi¬ 
tions  of  0.5  per  cent  H+,  0.5  per  cent  He+,  99  per  cent  Ot;  and  5  per  cent  H+,  95  per  cent  0+; 
for  temperatures  ranging  between  S00  and  1600°K, 


Fig.  6.  Dependence  of  ion-density  distribution  on  reference-level  composition;  T  =-1400°K. 
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Fig.  7.  Comparison  of  Blue  Scout  ascent  and  descent  data  with  theoretical  distributions 
calculated  for  T  =  1050°K  and  varying  reference-level  composition. 


n,{z) 


n,'(zo) 


exp 


go  Wi  E2(g  -  z0) 

*r  (E  +  z)(E  +  z„)_ 


Z  ».(*>)/ Z  ««(*<>) 

0  Q 


g0mq  R2(z  —  g0) 
fcT  (B  +  *)(B  +  zb) 


(7) 


Using  the  expression  lor  n,(z0)  given  by- 
equation  6,  the  distribution  of  the  ion  species 
can  also  be  written  in  terms  of  proportions  of 
all  constituents  at  the  reference  height. 


m,(z)  =  n,(z0)  exp 


go  wi,'E2(z  —  zp) 
kT(R  +  z)(R  +  z0)_ 


Z  a«(*o)  exp 


gow,  B2(z  —  z0)  1/2 

kT  (R  +  *)(B  +  z0)_ 

(70 


The  dependence  of  the  distribution  of  the 
ionic  species  given  by  equation  7'  on  tempera¬ 
ture  is  shown  in  Figure  8  for  a  reference-level 
composition  of  0.5  per  cent  IT,  0.5  per  cent  He+, 
99  per  cent  0*,  and  temperatures  varying  be¬ 
tween  800  and  1600°K.  These  curves  show,  as 
was  first  pointed  out  by  Mange  [1960,  1961], 


that  the  distributions  of  the  light  ions  will  reach 
a  maximum  at  some  altitude  above  the  refer¬ 
ence  level.  It  is  seen  that  increasing  tempera¬ 
ture  increases  the  altitude  of  maximum  density 
for  each  constituent.  Above  the  altitude  of  maxi¬ 
mum  density  an  increase  in  temperature  de¬ 
creases  the  rate  of  change  of  concentration  of 
the  ith  ion  with  altitude;  this  relationship  is 
reversed  below  the  peak. 

The  effect  of  varying  composition  is  seen  in 
Figure  9o  for  the  case  of  a  fixed  temperature 
of  1400°K.  In  Figure  96  normalized  distribu¬ 
tions  are  given  for  IT,  Ile+,  and  0+,  assuming 
the  same  proportions  as  in  9a  and  T  —  S00°K. 
For  fixed  temperature  and  reference-level  com¬ 
position  the  altitude  of  peak  density  for  each 
constituent  increases  with  decreasing  ion  mass. 
Also,  as  shown  in  Figures  9a  and  96,  the  alti¬ 
tude  of  maximum  density  for  a  given  ionic 
species  increases  with  an  increase  in  the  average 
ion  mass  at  the  reference  height.  The  curves  of 
these  figures  have  been  given  in  terms  of  nt  (z)  / 
nt(z0).  The  absolute  distribution  of  n,(z )  can 
be  obtained  by  multiplying  results  of  the  type 
given  in  the  figures  by  n.(z„),  the  experiment¬ 
ally  determined  density  of  the  ith  constituent 
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.  Fig.  S.  Dependence  of  the  vertical  profiles  of  the  ionic  constituents  O*,  He*,  and  H*  on  tem¬ 
perature  under  diffusive  equilibrium  conditions. 


at  the  reference  level.  The  shape  of  the  distribu¬ 
tions  given  in  Figures  9  will  not  change;  they 
will  merely  be  shifted  along  the  abscissa. 

As  indicated  by  equation  7'  and  by  the  results 
given  in  Figure  9,  the  scale  height  for  a  given 
species  is  a  complicated  function  of  emperature 
and  ion  composition  at  the  reference  level.  The 
procedure  outlined  above  provides  a  reliable 
method  of  determining  the  distribution  of  ions 


above  a  reference  level  under  conditions  where 
diffusive  equilibrium  exists.  The  proportion  of 
the  various  ionic  constituents  can  be  determined 
with  an  accuracy  of  about  15  per  cent  at  the 
reference  level.  An  upper  bound  can  be  given  to 
the  distribution  of  light  ions  with  a  relative 
abundance  of  2  per  cent  or  less  at  the  reference 
level.  The  solution  is  not  very  sensitive  to  varia¬ 
tions  in  composition  of  ions  of  mass  greater  than 


Fig.  9a.  Dependence  of  the  vertical  profiles  of  the  ionic  constituents  O*,  lie4,  and  H+  on 
reference-level  composition  for  T  =  1400°K. 
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Fig.  96.  Dependence  of  the  vertical  profiles  of  the  ionic  constituents  0+,  He*,  and  H*  on 
reference-level  composition  for  T  =  S00°K. 


16  amu,  since  the  relative  concentrations  of 
heavier  species  are  already  small  and  continue 
to  decrease  with  altitude. 

The  distributions  of  hydrogen,  oxygen,  and 
helium  ions  obtained  from  equation  7'  upon 
substituting  the  temperature  and  composition 
values  derived  above  for  the  Blue  Scout  flight 
?re  shown  in  Figure  10;  maximum  values  for 


the  helium-ion  distribution  deduced  from  these 
experimental  data  are  shown  as  curve  6;  He+ 
curve  a  is  the  distribution  obtained  for  1  per 
cent  He*  at  the  reference  level. 

These  results  show  that  H+  becomes  the  domi¬ 
nant  ion  at  650  km  and  that  the  maximum  in 
its  distribution  occurs  at  750  km.  The  helium- 
ion  concentration  reaches  a  maximum  at  625 


Fig.  10.  Distribution  of  O*  and  H*  with  altitude  and  the  upper  limit  for  the  He*  distribution 
obtained  from  the  Blue  Scout  results. 
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km,  200  km  above  the  reference  level,  and  is 
never  the  dominant  ion.  This  result  differs  con¬ 
siderably  from  the  conclusions  of  Hanson  [1962] 
and  Bourdeau  et  al.  [1962]  from  the  analysis  of 
daytime  measurements.  They  are  in  good  agree¬ 
ment,  however,  with  the  results  of  other  night¬ 
time  measurements.  On  the  same  nighttime  Blue 
Scout  flight  discussed  here  Ulwick  and  Pfister 
(private  communication,  1963)  measured  elec¬ 
tron  densities  by  means  of  an  impedance  probe 
experiment,  the  results  of  which  are  in  good 
agreement  with  our  positive  ion  density  meas¬ 
urements.  The  only  significant  difference  in  the 
data  from  the  two  experiments  occurs  below 
300  km  on  the  descent.  It  is  due  to  the  fact 
that  the  minimum  sensitivity  of  their  experi¬ 
ment  was  approximately  1.5  x  10‘/cm*.  The 
profile  of  electron  densities  given  by  Ulwick  and 
Pfister  [1963,  Figure  9]  is  a  composite  of  the 
ascent  and  descent  impedance  probe  data  ob¬ 
tained  on  the  Blue  Scout  flight  and  p-f  records. 
Donley  [1963]  reported  results  of  nighttime 
positive  ion  density  measurements  to  2700  km 
obtained  with  a  planar  electrostatic  probe  which 
are  also  in  good  agreement  with  the  results  pre¬ 
sented  here.  Preliminary  analysis  of  ion  mass 
spectrums  from  the  Ariel  satellite  [Willmore 
et  al.,  1963]  shows  large  differences  from  night 
to  day  in  the  mass  of  the  dominant  ion  above 
the  peak  of  the  F  region  consistent  with  the  re¬ 
sults  of  the  analysis  of  day  and  nighttime  ion- 
and  electron-density  profiles. 

The  experimental  results  and  the  analysis  in¬ 
dicate  that  helium  ions  must  be  present  in 
greatly  varying  amounts  in  the  lower  exo¬ 
spheres,  the  distribution  being  highly  dependent 
on  the  temperature  and  relative  concentrations 
of  all  the  ions  present.  This  is  in  agreement  with 
the  conclusions  of  Bauer  [1963]  and  Bourdeau 
and  Bauer  [1963]  as  to  the  temperature  de¬ 
pendence  of  the  helium-ion  distribution.  It 
should  be  noted,  however,  that  their  analysis  is 
based  on  a  special  solution  of  the  diffusive  equi¬ 
librium  equation  obtained  by  assuming  that 
there  is  a  transition  level  where  the  concentra¬ 
tions  of  two  ionic  species  are  equal.  This  as¬ 
sumption  is  not  implied  by  the  general  equation 
(equation  5). 

Taylor  et  al.  [1963],  for  example,  measured 
a  helium-ion  peak  at  about  600  km  and  con¬ 
cluded  on  the  basis  of  Bauer’s  analysis  that 
diffusive  equilibrium  did  not  exist,  since  the 


maximum  occurred  below  the  predicted  transi¬ 
tion  level.  Using  reference-level  values  obtained 
from  Taylor’s  data  at  450  km,  the  general  the¬ 
ory  predicts  a  broad  helium  distribution  with  a 
maximum  at  about  850  km  in  reasonable  agree¬ 
ment  with  Taylor's  experimental  results. 

( d )  Ion-Atom  Interchange  Coefficients 

The  determination  of  the  rate  coefficients  for 
the  processes  by  which  ionized  atomic  oxygen 
may  be  destroyed  is  important  to  the  under¬ 
standing  of  the  decay  of  the  nighttime  F  region. 
At  present  there  is  a  significant  disagreement 
between  the  values  deduced  from  theory,  lab¬ 
oratory  experiments,  and  ionospheric  data  for 
the  ion-atom  interchange  reactions  considered 
important  to  the  loss  of  0+  at  night  [ Bates  and 
Nicolet,  1960] : 


o+  +  n2  -» 

NO+  +  N 

(a) 

o+  +  o2  -» 

02+  +  O 

(6) 

The  theoretically  determined  rate  coefficients 
for  reactions  a  and  b  are  of  the  order  of  10'* 
cm*  sec'1  [ Field  et  al.,  1957 ;  Gioumousis  and 
Stevenson,  1958].  Lang stroth  and  Hasted  [1962] 
obtained  a  laboratory  value  of  the  rate  coeffi¬ 
cient  k,  for  reaction  a  above  of  4.7  X  10'“  cm* 
sec"1.  Sayers  and  Smith  [1963]  reported  a  lab¬ 
oratory  value  of  2.7  X  10"11  cm*  sec"1.  Talrose 
et  al.  [1962]  gave  an  upper  limit  to  k,  of  6.7  X 
10"“  cm*  sec'1.  Galli  et  al.  [1963]  published  a 
value  of  h,  =  2.2  X  10"u.  Using  the  Johnson  et 
al.  [1958]  ion  composition  data  obtained  from 
rockets,  Bates  and  Nicolet  [1960]  deduced  a 
value  of  kt  4-  0.1  fe,  ss  =  1.3  X  10"1*  cm* 
sec'1.  Yonezawa  et  al.  [1959]  obtained  a  value  of 
k,  =  3  X  10'“  cm*  sec'1  from  upper  atmosphere 
data  assuming  steady-state  conditions. 

Using  both  the  Blue  Scout  ion-density  profiles 
and  local  electron-density  distributions  provided 
by  G.  Stonehocker  of  the  National  Bureau  of 
Standards,  Boulder,  Colorado,  the  ion-atom  in¬ 
terchange  coefficient  k,  for  process  o  has  been 
evaluated  at  280  and  300  km  from  data  ob¬ 
tained  between  2030  and  0300  A.M. 

It  was  assumed  that  at  night  in  this  altitude 
range  there  are  no  significant  sources  of  ioniza¬ 
tion.  The  atomic  oxygen  ion  concentration  may 
change  as  a  result  of  radiative  recombination,  as 
a  result  of  ion-atom  interchange  with  Na,  O,, 
and  NO,  and  by  diffusion.  Because  the  NO  con- 
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centration  is  a  few  orders  of  magnitude  lower 
than  the  molecular  nitrogen  and  oxygen  concen¬ 
trations,  it  is  assumed  that  the  loss  of  0+  by 
ion-atom  interchange  with  NO  can  be  neglected. 
In  the  altitude  region  280  to  300  km,  Istomin’s 
[1960]  results  show  that  0+  contributes  more 
than  95  per  cent  to  the  total  positive-ion  con¬ 
centration  ;  0*  is  therefore  assumed  equal  to  the 
electron  density,  n.  The  rate  of  change  of  0+ 
must  then  be  approximately  equal  to  the  rate  of 
change  of  the  electron  density  with  time.  This 
may  be  expressed 

dn/dt  —  an 2  +  fejNgn 


+  &202n  +  di  v  (n  V)  (8) 
Dividing  equation  8  by  Nan, 


h  and  Zs,  are  the  same  order  of  magnitude,  and, 
since  (VN,  is  approximately  6  X  10“’,  y  ss  lcx. 
A  summary  of  the  calculated  and  measured 
quantities  used  in  the  evaluation  of  y  is  given 
in  Table  3.  dn/dt  as  dO*/dt  was  evaluated  using 
the  Blue  Scout  data  and  local  ionosonde  elec¬ 
tron-density  profiles  for  the  night  of  April  11-12. 
The  values  given  in  Table  3,  columns  2  and  4, 
are  time  averages  for  the  periods  shown,  a  was 
taken  equal  to  1.3  X  10'”  cm’  sec'1  [ Bates 
et  al.,  1939] .  The  divergence  term  may  be  writ¬ 


ten  [ Ferraro ,  1946;  Martyn,  1955;  Dungey, 
1956;  Ratclifje  and  Weekes,  I960]: 


div  (n  V)  —  Da  sin2  I 


c fn  . _ 3_  dn 

dz  +  2/7  dz 


where 


(10) 


I  —  magnetic  dip  angle  at  position  of  meas¬ 
urement. 

£>„  =  4.6  X  10'21  T/N  cma  see"1,  the  am- 
bipolar  diffusion  coefficient  [Dalgarno, 
1961], 

T  —  temperature,  °K. 

N  —  the  neutral  particle  density  cm"1.  U .  S. 
Standard  Atmosphere  [1962]  densities 
were  corrected  for  local  time  and  solar 
flux  conditions. 

z  =  altitude,  km. 

II  =  neutral  scale  height. 

The  values  of  n  and  its  first  and  second  deriva¬ 
tives  were  obtained  from  the  ion-  and  electron- 
density  profiles.  It  should  be  noted  that  the 
calculations  show  that  diffusion  acts  as  a  source 
below  300  km  and  changes  sign  slightly  above 
300  km. 

The  values  of  Na  and  O-  were  based  on  the 
work  of  Harris  and  Priester  [1963]  for  a  value 
of  the  10.7-cm  solar  flux  of  S9  watts/meter*  cps. 
A  small  correction  to  their  values  was  provided 
by  K.  Champion  (private  communication,  1963). 
As  shown  in  Table  3,  column  7,  kx  is  found  to 
vary  between  2  and  4  x  10'”  cm*  sec1.  There 


TABLE  3.  Calculated  Values  of  Ion-Atom  Interchange  Coefficients  (column  7) 
NO+  Production  Rates  (column  8) 


1 

Altitude, 

km 

2 

n, 

no.  /cm5 

3 

N* 

no.  /cm3 

4 

dn/dt, 
cm"3  sec"1 

5 

an", 

1  cur5  sec**1 

6 

—div  (nr), 
cm-3  sec-1 

7 

7> 

cm3  sec”1 

8 

ft,  0+N2, 
cm-3  sec-1 

280 

3.9  X  10s 

9.9  X  107 

55 

Time:  2030-2100 

0.2 

34 

2.3  X  10"12 

89 

300 

4.4  X  10s 

5.9  X  HP 

44 

0.25 

8 

2.0  X  10"12 

52 

280 

2.6  X  105 

8.4  X  HP 

16 

Time:  2100-0100 

0.1 

31 

2.1  X  10-» 

46 

300 

3.1  X  105 

4.8  X  107 

15 

0.1 

15 

2.0  X  10"11 

30 

280 

1.2  X  105 

6.8  X  107 

9 

Time:  0100-0300 
0.02 

24 

4.0  X  10"12 

32 

300 

1.9  X  105 

3.8  X  107 

S 

0.05 

20 

3.9  X  10“12 

28 
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is  an  uncertainty  of  about  a  factor  of  2  in  these 
calculations.  Nisbet  and  Quinn  [1963]  deduced 
a  value  from  the  loss  coefficient  /3  =  9  X  10'r' 
see'1  for  T  —  1050°K  at  300  km  from  a  statis¬ 
tical  analysis  of  ionosonde  data.  Substituting  a 
value  of  N,  =  4.8  X  107/cm3  (T  =  1050°) 
from  Table  3  leads  to  a  value  of  =  1.6  X  10'12 
cm5  sec'1.  Since  laboratory  experiments  in  gen¬ 
eral  provide  an  upper  limit  for  the  rate  coeffi¬ 
cients,  the  results  obtained  from  the  laboratory 
and  upper  atmospheric  data  are  in  reasonably 
good  agreement.  It  was  pointed  out  earlier  that 
nighttime  sources  of  ionization  have  been  neg¬ 
lected  in  these  calculations.  Hinteregger  (private 
communication,  1963)  has  estimated  that  the 
contribution  of  corpuscular  radiation  to  ion- 
pair  production  in  the  F  region  could  not  be 
greater  than  10  per  cent  of  the  daytime  solar 
ultraviolet  flux.  Including  this  ion  pair  produc¬ 
tion  rate  in  the  above  calculations  increases  the 
value  of  the  ion-atom  interchange  coefficient  k, 
by  a  factor  of  2. 

The  rate  of  formation  of  NO+  by  ion-atom 
interchange  process  a,  assuming  for  k,  the  val¬ 
ues  given  in  Table  3,  column  5,  gives  production 
rates  varying  between  28  and  89  per  cm3  sec'1 
(column  8) .  The  production  rate  decreases  with 
time  largely  because  of  time  variations  of  the 
O*  concentration.  This  represents  an  irreversi¬ 
ble  loss  of  nitrogen  from  the  upper  atmosphere. 
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